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Q ■ Abstract 

o ! 

The transfer matrix and matrix multiplication ansatz, when applied 
to nonequilibrium steady states in asymmetric exclusion processes and 
traffic models, has given many exact results for phase diagrams, bulk 
densities and fluxes, as well as density profiles and spatial and tempo- 
(T) ' ral correlation functions for models with a dynamics that is updated 

C^ i in (random) sequential and sublattice-parallel order. Here we consider 

fully parallel or synchronous dynamics, for which only partial results 
are known, due to the appearance of complicating strong short range 
qv ' correlations, that invalidate simple mean field approximations. 

This paper is based on two new ingredients: (i) a microscopic charac- 
terization of order parameters and local configurations in the relevant 
phases, based on the microdynamics of the model, and (ii) an im- 
proved mean field approximation, which neglects certain four point - 
and higher order correlation functions. It is conjectured that the den- 
q sity profiles, obtained here, are exact up to terms that are exponentially 

CJ , small in the system size. 
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1 Introduction 

Nonequilibrium stationary states (NESS) violate detailed balance, they can- 
not be described as Gibbs states, and their behavior shows a wealth of inter- 
esting phenomena that are absent in thermal equilibrium, such as boundary 
induced phase transitions, self organization, pattern formation, and long 
range spatial and temporal correlations. They occur in classical fluids |Q|, 
driven diffusive systems ^ ||, ||], granular flows || and lattice gas cellu- 
lar automata (LGCA) with collision rules violating detailed balance |0|, of 
which traffic models [j/J are simple examples. Unfortunately, a general the- 
ory like Gibbs statistical mechanics is lacking for NESS, where results seem 
much less universal, and depend strongly on boundary conditions, driving 
forces, and the (sequential or synchronous) order in which the microscopic 
dynamics is applied ||. 

The standard theoretical approaches are based on Langevin equations, 
fluctuating hydrodynamics, mode coupling theories, and ring kinetic the- 
ory, which are phenomenological and/or approximate in nature. A large 
amount of theoretical understanding has also been obtained from computer 
simulations as well. 

However since 1992 new methods for obtaining exact solutions for simple 
open one-dimensional systems, the so-called asymmetric exclusion processes, 
have been developed based on transfer matrix methods (Bethe ansatz, ma- 
trix product ansatz |9|- ||l9[| ). One can calculate bulk properties, phase dia- 
grams, density profiles of boundary layers and interfaces between coexisting 
bulk phases, as well as spatial and temporal correlation functions. The num- 
ber of exact results is rapidly growing. For an up to date list of references 
to analytic results for asymmetric exclusion processes we refer to the reader 
toRef. f§. 

Asymmetric exclusion processes can be characterized as open systems 
of interacting particles or random walkers on a lattice, in most cases linear 
chains, with hard core exclusion for double occupancy, and with hopping 
rates differing for different directions. The bulk dynamics may be deter- 
ministic or stochastic. Open systems are coupled to reservoirs at both ends 
through stochastic boundaries. Closed systems on a ring are coupled to lo- 
cal randomness such as defect sites [15, 21] or defect particles |l6|, [17], [H|. 



Applications range from shock waves in the Burgers equation, to traffic flow 
problems, reaction-diffusion systems and growth models. 

The problem can be formulated in terms of master equations with dis- 
crete or continuous time M , or in terms of the equivalent transfer matrices 



for spin chains |2C], or as microdynamic equations, as is usually done in 



LGCA, when used as models for nonequilibrium fluids [23]. LGCA rep- 
resent perhaps also a more faithful representation of real traffic problems 
[0, ^7j. The LGCA approach will be followed in this paper. 

The dynamics of updating sites may be applied in (random) sequential 
order, typical for the master equation description, or in parallel, i.e. fully 
synchronous for all sites, typical for LGCA, or in any intermediate version 
with strictly sequential or with sublattice-parallel updating ]1T| , 113]. The 
different ways of updating are an essential part of the model. They af- 
fect the existence of different phases in the phase diagram, as well as the 
structure of the spatial and temporal correlations. For instance, the so- 
called maximum current phase is present in the totally asymmetric exclu- 
sion process (TASEP) with open boundaries, when updating is carried out 
in random sequential order || 1C], but is absent for other updating schemes 



11, [0|, 21, 22]. The spatial correlations are weakest for random sequential 
updating, intermediate for sequential and sublattice-parallel updatings, and 
strongest for parallel updating. 

As far as analytic approaches are concerned many exact results about 
bulk properties, spatial and temporal correlations, and profiles are known for 
random sequential and sublattice-parallel updatings 0, |j| [| 11 1. The fully 
parallel updating schemes of LGCA offer the largest difficulties, because the 
dynamics creates strong short range correlations, which invalidate simple 
mean field approaches. Only some bulk properties like density, flux, and 



the phase diagram have been obtained [g, |21, |22j], where 'in the bulk' means 
' outside the boundary layers'. 

The main goal of the present paper is an analytic derivation of the cor- 
relation functions and the profiles, conjectured in Ref. [^] for an open 
TASEP with parallel updating, where the particles move forward fully syn- 
chronously at every time step with probability p = 1 (deterministic bulk 
dynamics) if the site in front of them is empty, and where input and re- 
moval rates specify the stochastic boundary conditions. In Ref. [|| a much 
richer stochastic version of the same model with p < 1 has been considered, 
and mean field results, based on a matrix multiplication ansatz, have been 
obtained. It includes the deterministic TASEP, discussed in this paper, as a 
special case, and it also contains the maximal current phase, even in the case 
of fully synchronous updating. However, even with the highly sophisticated 
matrix multiplication ansatz, the profiles and correlation functions have not 
yet been calculated as the associated matrix algebra is quite complicated 
0. In case of sequential or sublattice updating the model of Ref. M has 



already been generalized to include backward jumps with probability q, and 
input and removal rates at both ends of the chain fl2|| . 

Our approach starts in section 2 from the microdynamic equations, which 
describe the time evolution of the set of occupation numbers {tj(£)} of the 
sites i = 1,2,- ■ ■ , L on the chain at time t, as a discrete continuity equation, 
in which all variables are Boolean variables, having only values or 1: 

n(t + l) - n{t) = ji^it) - ji{t). (l) 

Here the instantaneous microscopic flux ji(t) counts the number of particles 
passing through link (i, i + 1) at time t = 0, 1, 2, • • •. The influx jo(t) and 
outflux J lif) specify the couplings to the stochastic reservoirs or blockage 
sites. Further specification of the fluxes depends on the model, and will be 
given later. 

A detailed analysis of the microscopic evolution equation for sets of par- 
ticle clusters TjTj + i • • •, derived in section 3, allows us not only to calculate 
the short range correlation functions that are built up through the dynamics, 
but also to identify the different phases, as well as the microscopic struc- 
ture of boundary layers and the interface between different bulk phases. 
In section 4 the phase diagram, bulk densities and fluxes in the different 
phases are calculated. The phase diagram shows a phase transition from 
a free flow regime (a < (3), where all particles are moving at maximum 
speed, to a congested or jammed regime (a > /?), where the dynamics is 
controlled by start-stop waves. Here a is the input rate and (5 the removal 
rate. When both rates are equal, there are coexisting phases with a sharp 
interface (shock wave) between them. The interpretation of the interface as 
a shock wave forms the direct link with Burgers nonlinear diffusion equation 



M 



In section 5 an exact hierarchy for particle-cluster correlation functions is 
derived, to which we apply our improved mean field approximation (MFA). 
It assumes for the low density phase, that inside the interface between bulk 
phase and boundary layer, the higher order correlations between on the 
one hand a particle-hole pair, and on the other hand the particle cluster 
directly in front of it, can be neglected. This is equivalent to neglecting 
certain 4-point correlation functions. For the high density phase, similar 
results are obtained through particle-hole symmetry. The MFA reduces the 
hierarchy to a set of recursion relations for the correlation functions and 
density profiles, which are solved in this section. Knowledge of the density 
profile also enables us to calculate the average travel time of particles. The 



results are new, and are in excellent agreement with extensive computer 
simulations. The present MFA is similar in spirit, but not in specific details, 
to the improved MFA of Ref. Jffl that accounts for short range 2-point 
correlations, but neglects higher order ones, and that leads to the solution 
of the TASEP on a ring without a blockage. 

To test the validity of our MFA we have applied the method in the 
appendix to the same open TASEP, but now with sublattice-parallel updat- 
ing, for which the exact correlations functions have been calculated in Refs. 



[11, 13]. It appears that the results agree with the exact results, apart from 
terms that are exponentially small in the length of the chain. Moreover, it 
turns out that spatial correlations arising from parallel updating are quite 
different from those coming from sublattice-parallel updating. 

In section 6 we exploit the analogy that a blockage site on a ring has 
with on the one hand the entrance site and on the other hand with the exit 
site of the corresponding open system. In this way we recover the results 



of Ref. [21] for the phase diagram, fluxes and the bulk densities for the 
TASEP with a blockage on a ring with parallel updating. In addition we are 
able to construct the higher order correlation functions, density profiles and 
finite size corrections (rounding) of the j (p)~velation ('equation of state') at 
the high and low density ends of the coexisting phase region. Again, the 
results are new and in good agreement with computer simulations. We end 
in section 7 with some conclusions and suggestions. 

2 TASEP with synchronous dynamics 

2.1 Definitions 

We consider a totally asymmetric exclusion process with open boundaries as 
a lattice gas cellular automaton, with particles living on a one-dimensional 
lattice with sites labeled i = {1, • • • , L}. The configuration of particles at 
time t = 0, 1,2, ••• is described by the set of occupation numbers {tj(£)} 
with i = 1,2, ••• ,L, where Tj(i) = 1 — at(t) = 1 if the site i is occupied 
by a particle, and Tj(i) = if the site is empty, i.e. is occupied by a hole 
(crj(£) = 1). The dynamics is defined such that all particles, which have an 
empty nearest neighbor site in front of them at time t, simultaneously jump 
to that site at the next time step (i+1). If that site is occupied, the particle 
does not move. So, the dynamics or updating in the bulk of the system is 
deterministic and fully parallel. 

Next, boundary conditions are specified. We consider an open system, 



coupled to two stochastic reservoirs, one that injects particles with a prob- 
ability a (0 < a < 1) into site 1, provided it is empty, and one that removes 
particles from site L with probability (3 (0 < /3 < 1) provided site L is 
occupied. 

If the system is in configuration {tj = Tj(i)} at time t, then, according to 
the dynamic rules described above, the configuration {t[ = tj(£ + I)} is given 
by the microdynamic equation r[ = Tj + jj_i — jj with bulk and boundary 
fluxes given byPI 

ji = Tjo-j+i (i = l,---,L-l) 

Jo = T ai = ao\ 

JL = t l g l+ i = (3t l . (2) 

Here a = tq and (3 = ctl+i represent a set of independent random Boolean 
variables, which take the values {0, 1} with expectations (a) = a and ((5) = 
(3, and which are drawn at every time step from a uniform distribution. 
For later analysis it is convenient to transform the microdynamic equation 
(|I])-(|2|) to hole-occupation numbers <Tj = 1 — Tj, yielding 

°"i = Tm+l + Vi-lCTi (i = 1,- ■■ ,L). (3) 

The equations for the time evolution of averages (Tj(t)) = (tj)j, correlation 
functions (riT^t etc. can be derived by multiplying the equations in (|l]) for 
different sites and subsequently averaging over arbitrary initial configura- 
tions {tj(0)}. In this way one obtains an open hierarchy which couples the 
time changes of a correlations function to higher order ones. 

For large times (t — ► oo) the system will approach a non-equilibrium 
stationary state, which is the main focus of attention in this paper. Averages 
over the NESS are denoted by (•••). This state is expected to be unique, 
i.e. independent of the initial configuration {tj(0)}. Therefore in analytical 
considerations the initial state is always taken to be the empty state {tj(0) = 
0} for all i {i = 1, 2, • • • , L), which is also the most relevant initial state in 
physical applications, such as traffic problems. 



In the model of Ref. [pi the boundary fluxes are the same, but the bulk flux is 
generalized to ji = piTtai+i, where Pi with < pi >= p represents a set of independent 
Boolean variables, similar to a and /3. 



2.2 Symmetries. 

A quantity of paramount interest is the average flux < ji > through the link 

(»,» + !), 

(ji) = {rm+i}, (4) 

and the local flow velocity or average speed vi = (Jj)/(tj). In the NESS these 
averages are independent of time, and the average flux can be calculated from 
the continuity equation, combined with (||) using the relation (tJ) = (tj). 
This yields a constant site-independent flux through the system, 

< ji-i >=< ji >= {nai+i) = j 
< jo >= «(1 - (n)) = j 

< JL >= P(tl) = j (5) 

with i = 0,1, ■■■,L. The flux in the NESS is translationally invariant. 
Once the nearest neighbor correlations are known, the density profile can 
be calculated from 

(Ti)=j + (TiT i+ x). (6) 

The equations of motion exhibit particle-hole symmetry, and the duality 
transformation, 

n *-► cr L - i+ i (i = l,---,L) 

t = a ++ a L+1 = J3 (7) 

Ji *~^> jL-ii 

maps the microdynamic equation (|l|)-(g) into the equivalent representa- 
tion (H), i.e. the microdynamic equation is invariant under particle-hole 
exchange. Consequently, the average occupation numbers satisfy the sym- 
metry relations 

(n}(a,/3) = (o- L -i+i)(/3,a) 

= l-{TL-4+l)(P,a) (8) 

with i = 1, 2, • • • , L. As the flux maps ji ^ Jl-% under the duality transfor- 
mation, the average flux satisfies the symmetry relation ji(a, /3) = jl-i((3, a). 
However, the average flux in the NESS is constant for all sites, hence 

j(a,P)=j(P,a). (9) 

The particle-hole symmetry is a very powerful tool, as all properties for 
a > /3 (high density) can be obtained from those for a < (3 (low density). 



3 Dynamics and Structures 

3.1 Build up of dynamic correlations. 

In this section we show that a qualitative analysis of the microdynamic 
equation and the resulting instantaneous configurations leads already to the 
complete phase diagram, to an identification of the relevant order parame- 
ters, and to a qualitative characterization of the structure of the high and 
low density phase, as well as to the typical dynamics in the different phases. 
In the sequel we will consider the dynamics of clusters of particles and 
holes described in terms of the Boolean variables, 

Tki = T fc T£; + i ■■■n 

Ski = CfcCTfc+i • • • a t (10) 

with k < I, i.e. a cluster has at least two constituents. The time evolution 
equations for these objects are obtained by multiplying (||) from k to I, using 
ajTj = 0. The result is 

Tki = T^i+i + Tk-i<JkTk+i,i+i 

= (u + Tk-\(Tk)Tk+i,i+i (11) 

Ski = Sk-l,l + Sk-l,l-lTl(Tl+l 
= Sk-ll-l((Tl + Ti<Ji + i). 



Multiplication of both equations in ( pT|) gives then 

(T kn S n+lj i)' = (r fc + r fc „ 1 o- fc )r A . +lin+1 S , n ^_ 1 (o-/ + 7707+1) = 0. (12) 

where the relation T n a n = has been used. 

The implications of (|12| ) are quite interesting, as it states that a config- 
uration containing (• • • 1100 • • •) cannot be created. As all possible configu- 
rations have evolved from the empty initial state, configurations containing 
a cluster of particles tailing a cluster of holes do not exist in the NESS. 

Moreover, a configuration (• • • 110100 • • •) with a single particle-hole pair 
separating the two clusters cannot be created either. The reason is that 
only the nonexistent configuration (• • -71100? • • •) in the previous time step 
could have created the configuration under consideration. The question 
mark represents a "0" or a "1". 

Similarly a configuration (• • • 11(01) 00 • • •) with k (k = 1, 2, 3, . . .) inter- 
mediate hole-particle pairs does not exist, as it could only have been created 



from the configuration (• • ?1 (10) fc 0? • • •) = (• • •?ll(01) fc - 1 00? • • •)• It then 
follows by complete induction that none of the above configurations can 
exist in the NESS. 

Consequently, the possible configurations generated by the dynamics 
from the empty initial state do not contain any configurations with two 
or more empty sites in front of the last cluster of particles. So, the fully 
parallel dynamics of the present TASEP builds up very strong short range 
correlations in the NESS. 

From the observations about the build up of dynamic correlations, we 
arrive at some important conclusions about the structure of the NESS. Let 
ko label the position of the tail particle in the last particle clustery. The 
configurations in the interval [1, k$ — 1] consist of isolated particles, separated 
by an arbitrary number of holes. In these so called free flow configurations 

Ti-m = (i < k ), (13) 

i.e. there is a 'hard core repulsion' between particles on nearest neighbor 
sites. The instantaneous fraction of occupied sites (density) in this interval 
is therefore p{< ko) < 1/2. The configurations in the interval [fco>-^] consist 
of isolated holes, separated by an arbitrary number of particles. In these so 
called jammed configurations 

o-i-io-i = (i> k ), (14) 

i.e. there is a 'hard core repulsion' between holes on nearest neighbor sites. 
The instantaneous density in this interval p{> ko) > 1/2. 

Therefore, large systems with an overall density (fraction of occupied 
sites) p < 1/2 necessarily have only bulk configurations of the free flow 
type with a narrow boundary layer of jammed configurations near the exit 
site of (yet unknown) width \r « L, whereas for p > 1/2 bulk sites 
contain only jammed configurations with a narrow boundary layer of free 
flow configurations of width A^ near the entrance site. Which value of p 
occurs, depends on the injection and removal rates, a and f3, and will be 
determined in section 4.1. 

3.2 Instantaneous profiles. 

To study the density profile we consider the dynamics of creation of pairs 
and larger clusters. We first observe that the first particle pair can only be 

4 Note, however, that an instantaneous configuration of the whole system may not 
contain any cluster of particles (at low densities) or any clusters of holes (at high densities). 



created at the exit site and only if (3 < 1. This follows from the evolution 
equation ( |ll|) for T^ k+i which shows that creation of a new pair at (k, k + l) 
requires the existence of a pair at (k + 1, k + 2). Therefore, creation of the 
first pair Tl-ij, at £ + 1 is governed by 

(tl-itlY = (tl-i + r L _ 2 a L _i) r L (1 - J3), (15) 

where the term tl_i(£)t£,(£) on the left hand side of the equation vanishes, as 
this pair has not yet been created. This implies that pairs cannot be created 
if {(3) = (3 = 1. Consequently, there is no boundary layer near the exit, the 
density profile is totally flat over the whole system and all configurations are 
pure free flow configurations. All space and time dependent correlations in 
this NESS can be calculated exactly in a simple manner, as will be shown 
in section 5.2. 

Next we consider the case /3 < 1, where particle clusters can be created, 
at least near the exit site. In case the injection rate is smaller than the 
removal rate, the average interval between arrivals 1/a at the pile up region 
near the exit is larger than the average interval 1/(3 between removals, and 
a large fraction of configurations are pure free flow configurations without 
any clusters near the exit site. On average there is only a narrow boundary 
layer of jammed configurations. So, the bulk properties of the system in 
the low density case a < (3 are determined by the injection rate a at the 
entrance site. 

In the jammed phase (a > (3 ), there is on average a large backup 
starting near the exit. The jammed configurations in the NESS cover the 
bulk of the system, leaving only a narrow boundary layer with free flow 
configurations near the entrance. The bulk properties in the high density 
phase are determined by the removal rate (3 at the exit site. 

In case a = (3, there exist downstream free flow configurations and up- 
stream jammed configurations, which occupy on average a finite fraction 
of the system. There are in fact two coexisting phases, spatially separated 
by an 'external driving field', imposed by the reservoirs with injection and 
removal rates a = (3. The phases are separated by an interface of micro- 
scopic size, containing only configurations with alternating particles and 
holes. These configurations, when present in a finite fraction of the system, 
would form the so called 'maximal current phase' ||. In the present model 
the width of this interface diverges, i.e. attains a macroscopic width, only 
at a single point in the phase diagram, where a = [3 — ► 1 (see section 4.2). 
In section 4.2 we will return to the dynamics and the structure of the in- 
terface. In the phase diagram of the stochastic model of [|| there is a line 
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(1 — a)(l — 0) = 1—p, where the density profile is flat over the whole system. 
In the present deterministic model (p = 1) this corresponds to the line (3 = 1 
(free flow phase), and to the line a = 1 (jammed phase). 

4 Phase Diagram 

4.1 Free flow and jammed phases 

For large system size L (thermodynamic limit) and a < (3 the system is in 
the free flow phase , and the dynamics rigorously implies that the bulk of the 
system (except the boundary layer near the exit) has only isolated particles. 
So, this low density phase is characterized by an average density pp < 1/2 
and has the vanishing order parameters (tjTj + i) = (TiTj + iTj + 2) = • • • = for 
i in the bulk phase, defined as i « L—\r with A# the width of the boundary 
layer near the exit. In fact, even the microscopic order parameter TjTj+i = 
in the free flow phase on the basis of (|13| ) . Of course the correlations (cr^cr^ +1 ) , 
(o"j(Tj_|_icrj + 2) etc. are non-vanishing in this phase. 

Moreover, the vanishing order parameter (tjTj + i) = for bulk sites, in 
combination with relations (||) imply the following properties of the low 
density free flow phase: 

J = in) = (r 2 ) = ••• = (n) = p 

(n) =3=9 = a/(l + «) 

(r L )=j/P = a/[P(l + a)], (16) 

where i is a bulk site. Several interesting features can be seen. First of 
all, the bulk dynamics is completely determined by the input rate a at the 
entrance site, as already explained in section 3.2. Secondly, the flux equals 
the bulk density, indicating that particles are never blocked in the free flow 
phase, and are traveling with an average speed vf = j/ < U >= 1. A 
particle entering the lattice is never blocked until it leaves the bulk and 
enters the boundary layer near the exit, where it slows down to a velocity 
V J = ]/{ t l) = (3 with (3 < 1 as a consequence of the pile up. The bulk 
density of the system, p = a/(l + a), is always smaller than 1/2, since 
a < (3 < 1. Therefore, we will also refer to the phase with a < [3 as the low 
density phase of the system. 

Similarly there is a jammed phase for a > (3 of high density p > 1/2, 
containing only isolated holes except in a boundary layer of width Ax, near 
the entrance. The phase is characterized by the microscopic order parameter 
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OiOi+\ = 0, or equivalently by the non- vanishing order parameters (tjTj+i), 
(tjTj + iTj + 2), etc. in the bulk. 

The properties of the high density phase with a > (3 can be related to 
those of the low density phase with a < (3 by the particle-hole symmetry of 
section 2.2. In this case, particle-hole symmetry implies that the dynamics of 
particles moving forward is identical to that of the holes moving backwards. 
In this point of view holes are injected at the exit site with probability (3 
and move downstream where they are finally removed from the lattice with 



probability a. So, (16) implies for the jammed phase: 

P = (T L )(a,/3) = {a 1 ){/3,a) 

= 1-<ti)O0, a) = 1/(1 + 0) (17) 

(tl) = (tz,_i) = ■■■ = (n) = p 

(n)(a,/3) = l-(T L )(J3,a) = l-P/[a(l + /?)], 

where i is a bulk site in the jammed phase, defined as i » Xl. Similarly 
the average flux in the jammed phase follows from (|9|) as 

3 = Jj(«, P) = 3f(P, a) = JLp = /3p=l-p, (18) 

which implies an average speed vj = j/p = (3. In the boundary layer near 
the entrance a particle has a higher average speed, v\ = j/(tx), and it slows 
down to speed /3, because of frequent blockage by a preceding particle. The 
average flux may also be written as j = 1 — p, indicating that the flux 
in the jammed phase equals the density of holes, which move with unit 
speed in backward direction. We note that these results are exact in the 
thermodynamic limit as L — > oo. 

In summary, the phase diagram, where < {a,f3} < 1, has a high 
density phase (p > 1/2) for a > (3, and a low density phase (p < 1/2) for 
a < (3. On the line a = (3 there are two coexisting phases, separated by 
an interface with an instantaneous position R: a low density phase can be 
found downstream of the interface and a high density phase upstream, each 
occupying a finite fraction of the system (as already discussed in section 3) 
with a jump Ap in the density between the two regions. In open systems 
with fixed a and (3 the position R wildly fluctuates, as it can be anywhere 
on the lattice with uniform probability. For a system of L = 1000 sites this 
statement has been verified by collecting 1.2 x 10 7 measured values of R into 
10 equal-sized bins. The resulting histogram is flat within fluctuations of 
1%. 



12 



Let us compare the behavior of flux j and bulk density p when crossing 
the transition line a = (3, where p = a/(l + a) and j = p in the low density 
phase, and p = 1/(1 + (3) and j = 1 — p in the high density phase. This 
shows that the flux j is continuous across the line a = (3, whereas the density 
makes a jump Ap = pj — pp = (1 — a)/(l + a). Therefore the NESS of this 
model shows a first order phase transition across the line a = [3. Figure 1 
shows the flux j(p), as a function of the average density at a fixed value of 
the removal rate (3. It is the analog of the equation of state for the pressure 
p(p) in thermal equilibrium. As a increases to /3, the flux j(p) = p, as well 
as the density increase up till pf = (3/(1+ {3). In the coexisting phase region 
(a = (3) the flux remains constant, and the bulk density is given by 

1 ^ (R\ a ( R\ 1 

It may be anywhere between pp and pj = 1/(1 + (3), as the location R 
hops around the interval (1,L). As a increases further ((3 < a < 1), the 
density increases from pj to 1, and the flux j(p) = 1 — p decreases. The 
triangle bounding the region (J < min{p, 1 — p}) is called the fundamental 
diagram in traffic problems, and corresponds to the coexistence region in 
thermodynamic phase transitions. As the interface position R in the open 
system, measured over a very long period of time, is uniformly distributed 
over all sites of the system, one obtains in the coexisting phase region a 
linear density profile p(x) by averaging the instantaneous profiles p(x\R) 
over all R, yielding 

fa) = \ E K*\B) = Pf{^ ~ I ) + PJj- (20) 

A more detailed discussion of the dynamics of the interface will be given in 
the next section. 

Before concluding this section we discuss the collective dynamics in the 
different bulk phases. In the low density or free flow phase every particle 
is preceded by at least one hole. Consequently it will advance one site 
per time step, and has the maximum speed vf = 1. In fact any free flow 
configuration with a hole in front of it, is propagated as a whole with velocity 
vp = 1. In the jammed phase on the other hand, every cluster of c particles 
is preceded by a single hole. At every time step only the lead particle of 
each cluster advances one site, and becomes the tail particle of the preceding 
cluster. Therefore clusters of constant length are moving backwards with 
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unit speed, and so do isolated holes. The c-th particle in a cluster makes its 
first move only after c time steps, and the average speed of a particle vj = j3 
(see (fj"8|)). The motion in the jammed phase is therefore characterized by 
start- stop waves, which are typical for most congested traffic flows. 

It is interesting to compare the results (|16|)-(|I8D for fully parallel updat- 
ing with those for sublattice-parallel updating, as derived in (A5)-(A7) of 
the Appendix. Let < r + > and < j + > be respectively the bulk density and 
the flux at even sites, as defined in the Appendix, and < j_ > and < r_ > 
the corresponding ones for odd sites, then we have for the sublattice-parallel 
dynamics in the free flow phase (a < (3), 

< t + > = a < r_ >= < t l >= a//3 (21) 

< j+ > = a < j_ >= 0, 

and for the jammed phase (a > (3) 

< r+ > = 1 < r_ >= 1 - <n>=l-(3/a (22) 

< 3+ > = P < J- >= o. 

The average properties, in particular of the odd sites, are very different from 
the TASEP with fully synchronous dynamics. 

4.2 Coexisting phases 

In this section we study the dynamics of the interface separating the low 
and high density phases in the coexisting phase region. To define the in- 
stantaneous location R and the instantaneous width w of the interface we 
introduce its front and tail sites 1$ and jo- Let jo be the position of the 
front hole in the first /most advanced cluster of holes, and Iq that of the tail 
particle in the last/least advanced cluster of particles, where jo < Iq (see 
section (3.1)). The sites jo and (.q belong with certainty to the low and high 
density phase respectively. Then the interval (jo, A) = jo + 2n + 1) contains 
only n alternating particle-hole pairs (10) n (n = 0, 1,2, • • •), which are al- 
lowed in both phases. The instantaneous position of the interface is defined 
as R = I (jo + 'o) an d its instantaneous width as w(n) = 2n = Iq — jo — 1. 
As the particle (c.q. hole) cluster is moving with unit speed forward (c.q. 
backward), the width decreases by 2 units per time step, vanishing after 
n times steps and yielding a hole-cluster adjacent to the particle-cluster. 
During this period R = ^(jo + ^o) remains fixed. 
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What happens next depends on the sizes h and c of the two adjacent 
hole- and particle clusters respectively. If h > c, the positions j'o and Iq 
remain fixed during (h — 2) time steps. After (h — 1) time steps the cluster 
of particles disappears, and j'o moves one site forward, whereas Iq makes 
a forward jump to the last particle on the then last cluster. If h < c, 
the same statements can be make with ^o/c/particles/forward interchanged 
with j'o//i/holes/backwards respectively. If c = h both points jo and Iq 
jump simultaneously. The dynamics of the interface width is illustrated in 
Figure 2. 

The probability distribution for the sizes of the right and left jumps, as 
well as those for the time intervals between the jumps, are determined by the 
probability distribution of finding n particle-hole pairs between two particle 
clusters in the high density region, or between two hole clusters in the low 
density region. The position R of the interface performs a random walk 
over all sites of the system around the average position (R) = L/2. Once 
the above probability has been calculated, its mean square displacement 
(SR) 2 = ((R — L/2) 2 ) and the associated short time diffusion coefficient T> 
can be calculated in principle for time intervals T satisfying the inequality 
5R = \J2VT « L/2. The long time diffusion coefficients vanishes due to 
the presence of the boundaries. Here we only illustrate the basic idea of the 
method by calculating the average width (w) of the interface, using a crude 
mean field estimate. 

The width (w), measured over a long time interval, is shown in Figure 
3, as a function of the injection rate. This behavior can be understood on 
the basis of simple mean field arguments. Let the instantaneous interface 
configuration be (• • •?00(10)™11? • • •) with (n = 0, 1, 2, • • •), then its width is 
w{n) = 2n. The probability on the configuration (00(10)™) tailing (11? • • •) 
- which is the start of the jammed phase - is P(n) = (1 — a)a n , where a n 
is the probability for injecting n particles and (1 — a) the probability for 
not injecting a particle. In the present asymmetric exclusion process, every 
injected particle is followed by a hole. The average interface width is then 

( w ) =^2nP(n) = 2/(1 -a). (23) 

n 

This estimate gives a fair estimate of the simulation results, as shown in 
Figure 3. 

5 Correlation Functions 
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5.1 Profiles and nearest neighbor correlations 

Consider first the low density or free flow phase where a < f3. To determine 
the density profile < r k > from (^) we need the nearest neighbor correlation 
function < T k r k+ \ >, which represents the probability of finding the sites 
(i, i+1) occupied. As already derived in section 4.1, we have < T k >= p and 

< T k r k+ i >= for bulk sites (k <C L — Xr). It remains to calculate these 
quantities in the boundary layer k ~ L — \r near the exit. To do so, we 
need the dynamics ([ll]) of the cluster correlation functions (10), averaged 
over the NESS, where < T" >=< T >. This yields 

< T ke >=< T k j +1 > + < Tfc_i<TfeT fc+1 ^ +1 >, (24) 

where I = k + 1, . . . , L and we recall that tl+i = I — [3- If T k+ \^ + \ refers 
to the least advanced particle cluster, then the sites (k — l,k) belong by 
definition to the boundary layer, and the probability for the configuration 
(tctT) k-i,i+i equals a times the probability for the configuration T k+ i^ + \. 
Since k-\-l is by definition the last site on the least advanced particle cluster, 
the occupation number a k equals unity with probability 1. In fact there are 
only very few particle clusters in the pile up region near the exit site, as the 
average removal interval 1/(3 is less than the average arrival interval 1/q at 
the pile up region. So, we expect that the least advanced particle cluster 
gives the dominant contribution to the probabilities, and we make the mean 
field assumption that the above factorization holds for all further advanced 
particle clusters as well, i.e. 

< T kl >=< T kHl >+a< T k+lil+1 > . (25) 

The present mean field approximation therefore assumes that 4-point corre- 
lations and higher order ones between a particle-hole pair and the particle 
cluster just in front of it (which by definition belong to the boundary layer 
in the present model) are negligible in the low density phase. 

The recursion relation above can be solved starting from / = L, where 

< T k) L + i >= < T k i > (1 — /?), and yields after iteration 

a(l-/3) <T f a(l-(3) \ L - k 
<T kL >= < T k+1:L > = I I <t l >, (26) 

where we have used the relation < T^l >=< tl > ■ Taking i = L— 1, L — 2, 
etc. gives < T k ^^i >= (a/ (3) < T k+ i j2 >, etc. and one finds by complete 
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induction 



a*' 



< T fc>L _ £ >= ( - ) < T k+l , L > . (27) 



Combining this with (^g) yields 



n \ L - k 



T U =< r k T k+1 . . . n >= [ - ) (1 - py- K <t l >. (28) 



The density profile is then obtained by inserting ( p8| ) for I = k + 1 in (n) 
with the result 

<^>_ i{ l + i^(|)" } , (29) 

where < r^ >= j//3 = a/ [0(1 + a)] has been used. By taking the thermody- 



namic limit of (28) and ( p9[ ) one recovers the properties of section 3 and 4 for 
the bulk phase, i.e. the order parameters < TkTk+i >=< TkTk+iTk+2, ■■ ■ > 
are vanishing, and the bulk density is p =< r^ >= j. In addition ( |29| ) 
describes the profile of the boundary layer near the exit. The excess den- 
sity, which equals the nearest neighbor correlation < T^Tk+i >, decreases 
exponentially on a length scale 1/ In ((3 /a), independent of the system size. 
Our mean field assumption, formulated below (p4[), also implies that site 
k in T^i marks the beginning of the boundary layer. Hence, the probability 
P(k) that the width of the boundary layer is L — k, is then proportional to 



the excess density < n > —j in (28), yielding after proper normalization, 

P(k) = c L - k /(i - C) (30) 

with 

C = a//3. (31) 

The expected width of the boundary layer in the low density phase is then, 
\ B = < L-k>=-±- = -^, (32) 

1 — C P — OL 



and the fluctuation 5Xb around this average is 



(8\ B ) 2 =< (L - kf > - <L-k > 2 = M 2__ = _ H (33) 



C _ a/3 
(1^"(M) 5 



The width of the boundary layer is 'microscopic' in nature, i.e. independent 
of the system size. In the limiting case as a j (3 (coexisting phases) the 
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width Xb diverges and becomes of macroscopic size. It is given by L — R, 
as the location R hops around over all sites of the lattice. 

The correlation functions in the high density or jammed phase with 
a > (3 can be obtained from particle-hole symmetry, e.g. < T k T k +\ > 
(a, 13) =< (JL-kVL-k+i > ((3, a) and yields 

< TkTk+1 >= TTp - TTp UJ • (34) 

Similarly, the density profile follows from (|34|) and @ as 



— >(£)' 



<^>=TT73[ 1 -< 1 - Q H«JJ' (35) 

It shows a boundary layer near the entrance site with a deficit density that 
decays on a length scale l/ln(a/ (3). 

For the coexisting phase region, where a = (3, the boundary layers at 
both ends diverge, i.e. become of macroscopic size, and the typical width £ 
of the interface £ =< w >= 2/(1 — a), or equivalently ~ 1/| lna|, as was 
calculated in section 4 on the basis of dynamical considerations. 

In Figure 4 the theoretical profiles for 2- and 3-point correlations are 
compared with those measured in computer simulations, for two different 
combinations of parameter sets {a = 0.50; {3 = 1} and {a = 0.50; (3 = 
0.51} and the results are indistinguishable. For (3 = 1 the profile is flat, as 



already explained below (15). Moreover we have compared simulations and 
theoretical predictions for the three- and four-point correlation functions 
< TfcTfc +1 rfc + 2 . . . > for a = 0.50 and j3 = 0.55. The former ones are shown 
in Figure 4b, and the results are again indistinguishable. 

In the Appendix the density profile and correlation functions for the 



same TASEP with sublattice-parallel updating [11] have been calculated, 
using the same mean field theory. Here the average flux is quite different, 
i.e. j = (3 < tl >= a, but the profile < r^ > /j is quite similar, i.e. 

f 1-/3 (a\ L ~ k ) 

< n >= j < <5fc,even H ~T~ I -5 J > , (36) 

at least for even sites. For comparison we also quote the results for the 

correlation functions, as derived in (All) of the appendix, i.e., 

a a 

<T kl > = <T k ^ 1 >=-<T k+1 , l >=-<T k+lil _ 1 > 



(1-P)?V-*) <TL >, (37) 

(3; 
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which should be compared with (|28|). We emphasize that k,l and L are 
even numbers in the present formulas. Comparison with (|2q ) shows that 
the correlation functions are quite different, in particular the differences 
between even and odd sites. The results (J36J) and (|37|) are exact to terms 
of order (a/ j3) L , as can be verified by comparison with the exact results 
derived in Refs. fll], 13 1. 



We conjecture that the corresponding results (Pq) - (|3q ) for the same 
TASEP with fully parallel updating and with open boundaries are exact as 
well, up to exponential terms of order (a/ ' f3) L for a < (3, and up to {(5/a) L 
for a > (3. However, for a close to (3, say a = (3(1 — 8/L), the expressions 
for the profiles in the boundary layers break down, as the neglected terms 
(a/ (3) L ~ exp(— 5) become of 0(1) for large L. 

5.2 Spatial and Temporal Correlations 

Next we will determine the correlation functions < TiTi+R > for bulk sites in 
the low-density phase (a < (3 < 1). As explained in section 4.1 the dynamics 
rigorously implies that the 'microscopic' order parameter TjTj + i = in the 
bulk phase. Consequently, the microdynamic equation for all sites outside 
the pile up region near the exit becomes 

t[ = r a 1 = d(l - Ti) 

t[ = n-i (Ki«L-\ R ). (38) 

The last relation implies translational invariance of the correlation for the 
bulk phase, i.e. 

< nn+R > = < n-iT i+R -i >= . . . = 

= < nn+R >= a < t r > -a < titr > . (39) 



The last equality follows from (|38D for t\. Then ( |39|) implies that the pair 
distribution function, defined as g{R) =< tiT\ + r > / < t\ + r >, satisfies 
the recursion relation 

g(R + l) = a-ag(R). (40) 

It follows straightforwardly from (|4(j) that the generating function has the 
form: 

^f^(«)= ';*-:> ■ (4i) 



R=0 
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By extracting the coefficient of z from ( |4~l| ) we find the pair correlation 
function in the bulk of the low density phase, 

< TiT i+R > = < t z >< r i+R > [1 + C(R)} 

= <Tl > 2 [l-(-a) R - 1 }. (42) 

This is an asymptotically exact result for 1 < {i, i + R} <C £ — Xr. In case 
j3 = 1 there is no right boundary layer, and ( f42|) is exact for all interparticle 
distances R. The pair function oscillates around its uncorrelated value, it 
vanishes for R = 1, as it should, and is independent of /3, because the bulk 
phase is only determined by the injection rate a. Neighbor particles are 
strongly correlated. The spatial correlation function in the bulk C(R) = 
— (—a) decays exponentially with a correlation length £p = l/|lna|. 
This correlation length becomes of macroscopic size as a — > 1 and = 1, 
where the system approaches the critical maximal current phase 0, which 
is in the present model only realized at the parameter values a = (3 = 1. 

The result (^) for the bulk correlations is very different from those in 
Refs. [11, [DJ for sublattice-parallel dynamics, where the spatial correlations 
do not depend on the interparticle distance R. For instance, in that model 
one has in the bulk of the low density phase 

_ / < n > 2 (1 + a) (i, R even) 
<r ^+*>-j o (elsewhere). (43) 

So far we have discussed spatial correlation functions for bulk sites of the 
NESS. Similar results can be obtained for time and space-dependent corre- 
lation functions. It follows from the microdynamic equation ( pq) that in the 
low density phase 

n +R (t) = n+R-i(t - 1) = . • • = T i+R -t(Q), (44) 

and consequently we have the asymptotically exact result, 

< T i+R (t)Ti(0) > = < Ti+R-tTi > 

= <r 1 > 2 {l-(-a)^- 1 }, (45) 

valid for < {i,j} <C L — Xr with j = i + R — t and R ^ t. For the special 
case (5 = 1 the result ((4EJ) is exact for all sites and for all times. 

One may also consider correlation functions involving fluxes. In the low 
density phase (a < (3) we have 

< T i+ R(t)n(o) >=< ji+ R (t)n(o) >=< j i+ n(t)ji(p) >> (46) 
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because the macroscopic order parameter TjTj + i = 0. The result obtained 
for the low density phase can be extended easily to the high density phase 
using particle-hole symmetry and the relation j = Oi + \ on account of (||) 



and (14). This yields for bulk sites 



< TiT i+R > (a, (3) =< 1 - T L _i + i - TL-R~i+1 - TL-i+lTL-R-i+l > (A Oi) 

1 r {i - (-#**} =< n > 2 {l - (-/J)^ 1 } 



{1 + 0? 

< nji+R > (a, (3) =< na i+R+ i > (a, j3) 

= < TL-R-i{l - T L -i+l) > (/?, a) 

= j <Ti >{l-(-(3) R+1 } 

< JiJi+R > (a, (3) =< a i+1 a i+R+1 >= j 2 {1 - (-/0)* -1 }. (47) 

For the coexistence region (a = (3) extensive measurements were made. It 
appears that the properties such as density, flux, etc. for both regions are 
equal to the properties in their associated phases. It seems interesting to 
measure the correlations between site i in the low density region and site 
i + R in the high density region. However, due to the random nature of 
the interface location, the occupation numbers Tj and Ti +R appear to be 
uncor related, even if R < l/|lna|. 

5.3 Profile on coexistence line 

In section 3 the profile and interface between low and high density phase 
(a = (3) has been studied on the basis of purely dynamical considerations. 
A few additional results can be obtained from particle hole symmetry (^) 
and the constant flux relation (||). 

On the transition line a = (3, it follows from (||) that < Tj >= 1— 
< TL—i+i >■ The density in the middle of the lattice is therefore on average 
1/2. The average number of particles < N > on the lattice follows by 
summing the above relation over all sites, and yields 

L 

<N>=Y / <r l >= \L. (48) 

i=\ 

The flux j = a/(l + a), which is continuous across the transition line (see 
section 3.1), allows us to calculate the boundary values exactly 

<n>=a/(l + a) ; < r L >= 1/(1 + a) (4g) 

< ri r 2 >=0 ; < cr L cr L _i >= 0. 

21 



There are two coexisting phases, downstream of the interface the low density 
phase with pp = a/(l + a) and upstream the high density phase with 
pj = 1/(1 + a), separated by an interface located at R. This interface hops 
around over the whole lattice, where two instantaneous profiles (at t\ and 
£2) are shown in Figure 5 for a = (3 = 0.50. 

5.4 Travel times 

This section studies the travel time of the particles, which is by definition 
the number of time steps that have passed from the moment a particle enters 
the lattice until it leaves the lattice. Because particles cannot travel with a 
speed larger than unity, the actual travel time will always be larger or equal 
to the size of the lattice T > L. For the average travel time, the following 
definition is used: 



<T>= 



L 

E 

8=1 



< n >/j = <N >/j, 



(50) 



where < N > denotes the average number of particles on the lattice. The 
resulting expression implies that the flux is equal to the average (unknown) 
number of particles < N > divided by the average travel time < T >. 

The present sections will give analytical and simulation results for travel 
times in low- and high density phases, as well as travel times in the coex- 
istence region. In the low-density phase of a large system (L — ► 00), the 
speed in the bulk equals unity. Therefore we expect the average travel time 
to be approximately equal to the size of the lattice, < T >~ L. It can be 
calculated directly from (|50|) and the density profile in ( p?9| ) and yields in 
the low density phase (a < /?), 



<T> 



L + 



L + 



1-/9 
j3 — a 

1-/3 

/? 



a 



1 



(L 



00 



(51) 



where the term (a//3) L should be neglected for consistency. We observe that, 
for a and j3 not to close to one another, the travel times are of the order of 
L, as expected. Moreover, the average travel time < T > is a function of 
both parameters a and j3. As a j 0, the neglected correction terms (a//?) 
in (J5l|) and ( p9[ ) become of 0(1), and ( |5l| ) is no longer valid. 
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For a large range of combinations of a and (3, the average travel times 
have been measured with the simulation program, which is able to measure 
the travel times of a variable set of individual (labeled) particles. The mea- 
sured travel times are in excellent agreement with relation (|5l|), even very 
close to the transition line (a = (3). In Figure 6a we show a histogram of 
travel times in a system with L = 1000, a = 0.50, (3 = 0.51 in the low 
density phase near the transition line a = (3, where travel times may vary 
considerably. The average travel time for the specific measurement was 1049 
time steps. Expression (pj) gives < T >= 1049, so the agreement is very 
good. 

In the high density phase (a > (3) the average travel time should be 
approximately L//3, because the speed in the bulk of the system is equal 
to vj = (3. The average travel time < T > can again be determined from 
definition ( |50| ) and ([35]) with the result 

<T> =-R-^R (£-«>). (52) 

p a — (3 

For a = 1 the average travel time assumes the value L//3. For a < 1, the 
average travel time is smaller than L//3, because then the particles travel 
with a velocity larger than (3 just after entering the lattice. Comparison with 
the results of the low-density phase also shows that the statistical spread 
in travel times is larger in the high-density phase, because the particles 
are often blocked by other particles. Figure 6b shows a histogram of the 
measured travel times for a lattice of L = 1000 sites, and injection and 
removal rates of a = 0.51 and [3 = 0.50 respectively with a predicted average 
< T >= 1951. Obviously there is a large spread around the average travel 
time (|5l"|). The distribution is not symmetric and therefore not a Gaussian. 
The asymmetry of the distributions in Figures 6a, b are related through 
particle-hole symmetry. 

The average travel time < T > on the transition line a = (3 can be easily 
determined, using the average number of particles in (|50|). The travel time 
combined with the flux j = a/(l + a) yields then 

< T >= < N >/j = \L (l + I) . (53) 

This equation merely shows that in the coexistence region, the average travel 
time is just the average of the low- and high-density travel times, because 
both phases are present and occupy on average an equal fraction of the 
system. The long time average < T > agrees very well with (|53[) . Individual 

23 



particles will have a considerably shorter travel time of order L when R is 
located near the exit site, and a longer one of order L/a when R is located 
near the entrance site. 

The prediction ( p3[) is in excellent agreement with the simulation re- 
sults, as shown in Figure 7. We also note that the theoretical description 
of our totally asymmetric exclusion process in terms of occupation numbers 
(indistinguishable particles) does not allow us to calculate travel times of la- 
beled particles, or calculate the probability distributions in Figure 6. Such 
calculations would only be possible in a TASEP with labeled particles, as 
considered in Refs. [17, 19]. 



6 TASEP on a ring with blockage site 

The fully synchronous TASEP of section 2.1 as a closed system with N = pL 
particles and obeying periodic boundary conditions is a fully deterministic, 
rather uninteresting system in its NESS. For p < ^ all sites are in pure 
free-flow configurations and particles travel, say, counterclockwise with unit 
speed and flux j = p. For p > ^ all sites are in jammed configurations, and 
holes are traveling clockwise with unit speed, and the flux j = 1 — p. 

The dynamics becomes more interesting by inserting a stochastic block- 
age at site i = L with a transmission rate j3 < 1. The microdynamic 
equation for sites i = 2, 3, . . . , L — 1 is the same as in (|l|)-@, but the fluxes 
referring to the blockage sites are 

k =3L = /3t l cti, (54) 

where the Boolean variable (3 with expectation < (3 >= /3, is defined in 
a similar manner as a and j3 below (0). For j3 = 1 one recovers the fully 
deterministic case with periodic boundary conditions. 

First we observe that the dynamics at fixed (3 is invariant under the 
duality transformation 

Ti <-► (JL-i+l 

p<->l-p, (55) 

and that the average occupation satisfies 

< n > (p, 13) =< <TL-i+i > (1 - p, 13). (56) 

A mean field theory for the bulk properties of this model in the thermody- 
namic limit has already been given by Yukawa et al pi]], as well as extensive 
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numerical simulations, specially for the coexisting phase region. However, 
the correlation functions and density profiles have not yet been studied ana- 
lytically. In the Appendix we have discussed the TASEP with open bound- 
aries and sublattice-parallel updating. The corresponding models on a ring 
with a single blockage site and sublattice-parallel updating or with random 



sequential updating have also been solved exactly in Refs. |l4|-[|ljj]. How- 
ever, the analytic results show little similarity with those for the present 
model, and will not be discussed further. 

The build up of dynamic correlations and structures may be analyzed in 
a similar way as in section 3 for open systems, and one recovers (15), except 



for site i = L. This implies that a cluster of holes upstream of a particle 
cluster can only be created at the blockage site. Of course, hole clusters 
upstream of particle clusters may be present in the initial state anywhere on 
the lattice. In the low density phase, such configurations will be destroyed 
in a time that is roughly equal to twice the size of the largest hole cluster. 
Moreover, from a detailed analysis of the dynamics, similar to section 3.1, 
one can derive exact relations for the microscopic dynamic correlations for 
the bulk of the low density phase, such as 

TjTj + i = (bulksites i <^ L — Xb) 
t~lt~i = (blockage site), (57) 

where Xb is the width of the pile up region, downstream of the blockage. 

The continuity equation (jl]) and (]54|) yields then in combination with 
( |57| ) for the low density phase in the NESS, 

j = (3 <t l > = < n >=< r 2 >=...=< n > {i<^L- x B ) 

= < n > - < Tm+i > (L - X B ~ i < L). (58) 

In the low density phase there is an excess density in the pile up region. 
Consequently, as the total density p = N/L is fixed, the density at bulk 
sites has the form < t\ >~ p{l — 0(Xb/L)}, as we shall see later. 

In fact, one can infer most of the results for the ring model with a 
blockage from section 4, by considering the flux j = (3 < tl > across the 
link (L, 1) as the influx jo appearing in (||) for the open system. This relation 
defines the effective input rate a e through the relation j = a e (l— < n >) 
and yields in combination with j =< T\ > in (|5^ 



a e = j/(l-j) (59) 
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where a e approaches p/(l — p) in the thermodynamic limit. The flux across 
the link (L, 1) can be equally considered as the outflux < jl >= f3 e < tl> 
in (S) of the corresponding open system. This identifies the effective removal 
rate 3 e = 8 as the transmission coefficient of the blockage site L. 

The phase diagram for the system with a blockage can then be read 
off from Figure 1, showing the flux j(p) of the open system at a fixed re- 
moval rate 8. Consequently, for a e < /3 e , or equivalently for j < j3/{l + j3) 
(where j ~ p for large systems) , the system is in the low density or free flow 
phase. If the density p approaches pf = 8/(1 + (3), or if a e approaches 
/3, then the system enters the region of coexisting phases, and the pile 
up region, which had before a microscopic width A# of approximate size 
l/ln(/3 e /a e ) = l/ln[/3(l — p)/p], grows to macroscopic size, as in a wetting 
transition. For p > pp an interface appears downstream of the blockage, at 
a location R, and the pile up region has the macroscopic size L — R. As the 
density p increases further to pj = 1/(1 + /?), the location R moves further 
downstream, according to (see (|T^)) 

R\ 8 / R\ 1 

+ 1-t )-r-is (pf<p<pj). (60) 



LJ 1 + 3 V LJ l + B 

As p | pj, the free flow phase disappears (R — > 0), and the system goes into 

a pure jammed phase, where there is again a microscopic boundary layer, 

just upstream of the blockage with a deficit density. 

In summary, the TASEP on a ring with a blockage site has the following 

phases, 

• free flow phase: p < pf = /3/(l + 8) 

• coexisting phases: pf < p < p,j 

• jammed phase: p > pj = 1/(1 + 8). 

The above results were first obtained and verified against computer simula- 



tions by Yukawa et al. [21]. 



The density profile < Tj > in the free flow phase can be inferred from 
the corresponding profile (|29|) for a e < 8 e , and yields, 

<n> = j{l + ^C L ~ 1 } (i<L) 

<t l > = j/8, (61) 
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where terms of 0((, ) have been neglected. The relation above is valid for 

c= f-«w) <1 or 1<PF - <62) 

In fact, the first line in (|6l| ) also covers the case i = L. The relation between 
flux j(/3) =< n > and density /? in the free flow phase follows by summing 

where 0(£ )-terms have been neglected. The flux j(p) can be solved from 
this quadratic equation, where the root with the minus sign is the physical 
root. For large systems j differs only slightly from p. However, the (D(l/L)- 
correction becomes more and more important as j f pp = (3/(1 + (3) where 
the denominator in (p3|) diverges. By a perturbation expansion to 0(1/ L) 




we find, 
with 

1 - p\ 1 _ Ap 

.TTp) I = T' (65) 

where Ap = pj — pf is the difference in density between the two coexisting 
phases. The numerical solution of ( |63| ) is plotted in Figure 8 as the dashed 
line. 

We have again performed computer simulations on large and small sys- 
tems to test the density dependence of the flux j(p). After preparing the 
system in a random initial configuration, we let the system relax for 4 x 10 4 
time step, after which it is assumed to be in the NESS. We have calculated 
time averages over 3 x 10 5 time steps, and ensemble averages over 50 dif- 
ferent initializations. The agreement in the interval p < pf between theory 
and simulations is excellent, even for small L and 0, and for p close to pf 
where the difference between j and p is largest, as shown in Table 1 and 
Figure 8. For pp < p < \ the difference between theory and simulations 
becomes somewhat larger. The reason is that (qTJ)-(^) are only valid for 
( < 1 or j < pp. As soon as £ ~ 1 or j ~ pp equation (|63|) starts to loose 
its validity because in ( pl[) and (|||) terms of 0(( L ) ~ 0(1) have been ne- 
glected. Simulation results for j(p), similar to those in Figure 8, have been 



presented in Ref. [21] without a theoretical explanation. 
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Next we consider the n-point correlations, which are given through (|2 
and (lilf), i.e. 

< r k r k+1 . . . r k+n >= (i//3) (1 - (3) n C L ~\ (66) 

with £ given in (p3). This relation also gives the profile < Tf. >= j+ 
< TfeTfc + i > of the pile up region downstream of the blockage. In Figure 9 
the n-point correlations with n = 1,2 have been compared with computer 
simulations and again there is good agreement. 

Moreover, through arguments similar to those in section 5.1, we conclude 
that the probability P(k) to find the first site of the pile up region at site k is 
proportional to the excess density < T k T k +i >, so that P(k) = ( L ~ k /(1 — (). 
The average width A# of the blockage region and the fluctuation 5Xb around 
this average are then found from ( |32] ) and (|33|) with ( given through (p2|), 
i.e. 

X B = <(L-k)>=j/[(3-j(l + (3)} 
(5\ B ) 2 = < (L-k) 2 >-<L-k> 2 

= Pj(l-j)/lP-j(l + P)}\ (67) 

where < . . . >= J2k(- ■ -)P(k)- Simulation results for A# and 5\b have been 
presented in Ref. pi] ]. 

For the high density phase all corresponding results can be obtained from 
particle-hole symmetry. For instance, the profile of the depletion region 
just upstream of the blockage site L is given by (pS|) with a e = (3 and 
Pe = i/(l — J) where j(p, /?) in the jammed phase equals j(l — p, (3) in ( |63| ) 
in the free flow phase. 

The behavior of the interface in the coexistence region is very different 
from that in the TASEP with open boundaries. In the latter the location R 
of the interface for a given a = (3, can be anywhere on the lattice with equal 
probability, as the actual density fluctuates between (3/{l + j3) and 1/(1 + /?). 
In the TASEP on the ring with blockage the density p = N/L is fixed, and 
R is on average given through (^0|), and there are only small fluctuations 
5R around the average it!. In summary, the previous discussion confirms our 
intuitive interpretation that the TASEP on a ring with a blockage behaves 
the same as the TASEP with open boundaries with injection rate a e = 
j/(l — j) and removal rate j3 e = (3. This similarity holds for bulk properties, 
profiles and correlation functions. 
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7 Conclusion 

In this paper, which is in part an account of Ref. p2| , we have studied the 
non-equilibrium stationary state (NESS) of the totally asymmetric exclusion 
process (TASEP) with fully synchronous and deterministic (p = 1) bulk 
dynamics (i) for open systems, coupled to particle reservoirs with injection 
rate a and removal rate j3, and (ii) for closed systems on a ring containing 
a stochastic blockage site with transmission rate (3. Mean field theories 
and the Boltzmann equation give a totally inadequate description of these 
far-from-equilibrium states, because of the existence of strong short range 
correlations. 

The theory presented here is based on two new ideas: (i) we derive, 
starting from the microdynamic equations for the TASEP, the explicit mi- 
croscopic specifications of the configurations and order parameters for the 
separate phases; and (ii) we introduce an improved mean field approxima- 
tion in ( Pq ) that neglects fourth and higher order correlation functions at 
the interface between the bulk phase and the boundary layers. The results 
for the profiles have been compared with extensive computer simulations, 
and turn out to be indistinguishable from the analytic results. We therefore 
conjecture that our results for the open TASEP with a ^ (3 are exact up 
to terms that are exponentially small in the system size L, for instance of 
order (a/ [3) L for a < [3. Our results for the TASEP on the ring with a 
blockage show small differences between theory and simulations. Clearly 
the identification of the flux through a blockage as both the influx (to define 
a e ) and the outflux (to define j3 e ) of the open system is only approximate. 
Moreover, the neglected correction terms (a//3) L start to become of 0(1) as 
a e t Pe or j | Pf (see Figure 8). 

The first idea has enabled us to obtain exact results not only for bulk 
densities and currents, but also for the spatial and temporal correlation 
functions. The second idea has enabled us to obtain analytic results for the 
profiles in the boundary layers of density < r^ > and cluster correlation 
functions < r^rfc+i • • • >. For the more general stochastic model with p < 1 
of H no analytic results for profiles and correlation functions are known. 
The ideas in (i) are akin to the elimination of the 'Garden of Eden' states 
in ||, and those in (ii) to the 'paradisical mean field approximation', hinted 
at in Ref. ||, but that lingo is not ours. 

It is of interest to compare our results for fluxes and bulk densities with 
known results. The phase diagram has been obtained before in [11, [D|] for 



the same TASEP with sublattice-parallel updating, and in [21] for the fully 
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synchronous TASEP on a ring with a blockage. It has a free flow phase 
(a < (3), and a congested phase (a < (3), which coexist when a = j3. 

In Ref. pi a more general stochastic TASEP with fully synchronous 
dynamics has been analyzed, where particles hop only with probability p 
(with p < 1). The corresponding phase diagram contains the present {p = 1) 
phase diagram for a and (3 less than a c = 1 — yjl — p), but it also contains 
more phases, such as the maximal current phase. In our {p = l)-model the 
maximal current phase occurs only as a bulk phase for the special parameter 
values a = (3 = 1. In addition, the interface of finite width, separating 
the coexisting phases for a = [3 < 1, constitutes a 'microphase' of maximal 
current configurations. In the TASEP with p = 1 the interface region (see 
section 4.2) contains only local configurations common to both the low and 
high density phase. These common local configurations are identical to the 
local configurations that constitute the maximal current phase. It would be 
interesting to find out for the fully synchronous TASEP with p < 1 if it is 
possible to identify an interface region for a = (3 < a c , which constitutes a 
'microphase of maximal current configurations' as well. 

In the low density regime (a < [3 < 1 — ^/l — p) of Ref. |8| the flux and 
bulk density are found as 

p-a a(l-a) 

j = a 7T and p = ^-, (68) 

p — a z p — a z 

which reduce for p = 1 to the results in (16). The corresponding properties of 
the phase diagram in the high density regime can be obtained from particle- 
hole symmetry. 

In order to illustrate the particle-hole attraction in these models, and 
its dependence on the hopping rate p, we compare our results for the deter- 
ministic version (p = 1) for the nearest neighbor correlation functions in the 
bulk (for i <^ L — \r), i.e. 

< TiT i+1 > = ; < o-j<Ti + i >= 1 - 2p (ie bulk) 

< na i+1 > = < am+t >= j = p (Vi) (69) 

with those for the stochastic version (p < 1) in Ref. ||, reading 

< rjTj+1 > = p - j/p (i e bulk) 

< o-j<7j + i > = 1 - p - j/p (i e bulk) 

< Tj£Tj +1 > = <am + i>=j/p (Vi), (70) 
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where relation (6) has been used. Of course ( |70| ) includes ( p9| ) for p = 1. 

We notice that in the low density phase of the deterministic version 
there is a 'hard core repulsion' for nearest neighbor sites, as < TjTj + i >= 
or equivalently, a strong attraction of particle-hole pairs on nearest neighbor 
sites. In the stochastic version (p < 1) there is also a particle- hole attraction 
on nearest neighbor sites, because the covariance 

a 2 (p — a) 2 

< na i+1 > - <n>< a i+1 >= —± £- > 0, (71) 

p(P — a A ) 

i.e. there exists a positive correlation between an occupied site and the 
empty site, just in front of it, which increases monotonically as p j 1. 

It would be very valuable to extend the present method, based on the 
ideas summarized in (i) and (ii) at the start of this section, to calculate 
the profiles, and spatial and temporal correlation functions in the stochastic 
TASEP of Ref. || with its much richer phase diagram containing a maximal 
current phase. 

Appendix A: Sublattice— Parallel Dynamics 

In this appendix we will illustrate how the method of the present paper, when 
applied to the TASEP with sublattice - parallel updating, yields the exact 
results for the bulk properties and the asymptotic (large system) results for 



profiles and correlation functions, as obtained in Refs. [11, 13 1 The dynamics 
consists of two substeps. In the first step from t — ► t' = t + s, the pairs 
(1, L), (23), ... (L — 2, L — 1) are updated in parallel where L is even. There 
is only a possibly non-vanishing microscopic flux % through the even link 
(i, i + 1), whereas the flux ji through the odd link vanishes. In the second 
step, from t ->■ t" = t' + \ = t + 1, the pairs (1,2), (3,4), . . . (L - 1, L) are 
updated in parallel, then ji + \ is possibly non-vanishing and % vanishes. 

If we denote T a (t),T a (t + |) and r a (t + 1) respectively by T a ,r' a and r" 
with a = {1,2,... ,L}, then the microdynamic equation for the first step 
t — ► t' becomes: 



t[ = TiTi + \ (0 < i < L;i even) 

r i+i = T m + T i a i+i (0<i < L;i even) 

and for the second step t' — ► t + 1 

^i =< + <-!< (0<i<L;ieven) 

T i+i = T Ui T t+2 (0<i<L;i even). 
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(A.l) 



(A.2) 



With the conventions tq = a and ctl+x = P, as defined below (|2|), these 
equations also include the boundary conditions for the open system. More- 
over, we observe that the evolution equations (A1)-(A2) are invariant under 
particle-hole exchange. 

In the NESS there is again a constant site-independent flux through the 
system. The flux < ji > out of even sites at integer times is equal to the 
flux < j[ +l > out of odd sites at half-integer times, i.e. 

j = < Ji > = < TiOi+l > 

= < j' j+1 >=< (r i+ i + Tia i+1 ){a i+ 2 + T i+2 o- i+3 ) >, (A.3) 

where tq = a and o~l = (3 and i is even. 

By studying the dynamics of clusters as in section 3 one finds that (Al)- 
(A2) impose some very strong constraints on the allowed configurations, 
i.e. 

ira)'U,i = rUr'A-A = 

{tTO-o) "_ M+2 = T 'i-l T 'i (T 'i+l (T 'i+2 = 
{(TTT)i- 1)i+1 + (TaTT)i-2,i+l}0-i0-i + lT i+2 Ti +3 = . (A.4) 

This implies that configurations containing (. . . 10+ . . .) and (. . . 1100+ . . .) 
can never be created at integer times if one starts from an empty initial 
state. Such configurations are therefore absent in the NESS for any value 
of a and (5. The subscripts (±) on r± indicate that the relevant site has an 
even (+) or an odd (— ) label. Of course configurations (. . . 10_ . . .) as well 
as (. . . 0011 . . .) are allowed at integer times. 

Moreover, by arguments similar to those in section 3.1 and 3.2, one shows 
again that the first particle-cluster can only be created at the exit site. Let 
ko be the location of the last particle on the last cluster of particles, then 
the interval downstream of ko contains only isolated particles, separated by 
holes (free flow configurations with interval density p(< ko) < i), and that 
upstream of k^ contains only isolated holes (jammed configurations with 
interval density p(> ko) > o)- All conclusions in the last three paragraphs 
of section 3.1 carry over to the TASEP with sublattice - parallel dynamics, 
and so does the phase diagram. 

Next, we consider the bulk properties for large systems (L — ► oo) in 
the free flow phase (a < j3 and pf < V)- The low density phase is again 
characterized by the microscopic order parameter TjTj + i = 0. To perform 
these calculations, we start, as in section 4.1, from the constant flux relations 
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(A3) in combination with the vanishing order parameter, 

3 = < Ji >=< t 2 >=...=< n >= (3 < t l > 
= <3i+i >=a(l- <n >) =<t 2 + t 3 >= 
= ...=<Ti + Ti+i > (A.5) 

valid for i <C L — A# and i even. This implies for odd bulk sites < r_ >= 0, 
and for even bulk sites < r + >= a and < tl >= a/(3. If < j+ > denotes 
the flux out of even bulk sites, and < j'_ > the one out of odd sites, then 
the results for the free flow phase (a < (3) can be summarized as, 

< r + >= a < r_ >= < t l >= a//3 

< j+ >= a < j. >= . (A.6) 

The corresponding relations for the jammed phase (a > (3) can be obtained 
from the relations (B) for particle-hole symmetry, and read: 



(A.7) 



(A.8) 

Next, we consider the profile in the right boundary layer of the low density 
phase, and we construct the dynamics of the cluster functions T, as in flllD. 
By specializing these equations to the low density phase one arrives after 
lengthy, but straightforward algebra, at a coupled hierarchy of equations for 
the correlation functions. 

Let i or i + 1 (with i = even) be the last particle position on the last 
particle cluster, then we find in the NESS the exact relations, 

< T ik > = < T i)k _i >=< T ijfe+1 > + 

< (rcro-r) i _ 2 ,fe+i > + < {TO-TO-T)i_ 2 ,k+l > 
< T i+1}k > = < T i+1>k _x >=< T i+ljfe+ i > + < (rcrT) iifc+1 >, (A.9) 

where i and k are both even. The present set of coupled equations is the 
analog of (|24|). By application of the mean field assumption, formulated 
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< t + > (a, /3) = 1- 


- < t- > {j3, a) =1 


< t_ > (a,(3) = 1- 


- < t + > ((3, a) =1-/3 


< ri > {a, 13) = 1- 


- < t l > (/?, a) = 1 - (3 /a, 


and for the fluxes, using (^), 




< 3+ > (a, (3) 


= <3+> (13, a)=p 


<3- > (a, (3) 


= 0. 



in (|24|)-(p5|), the above set of equations simplifies to the set of recursion 
relations with i and k even, 

<T ik > = < T i>fc+ i > +a < li+i.fe+i > +a < T i+2i k+i > 
< T i+ljk > = < T i+1)k+1 > +a < T i+2 ,k+i > ■ (A.10) 

The boundary condition for this set is included by setting k = L. 

Our special mean field assumption for the low density phase neglects 
again higher order correlations between on the one hand the particle-hole 
pairs to" and to~to~, and on the other hand the tailing particle cluster at the 
interface of bulk phase and boundary layer. The solution of these recursion 
relations yields 

y i+ \i- p)U^). (A.n) 

From these results and from (y) we obtain the density profile of the boundary 
layer near the exit for even sites i, 

<Tl >=a + (l-P)^-j . (A.12) 

For odd sites i + 1 follows similarly, 

<r m >=(l-/?)(^) L l ■ (A.13) 

The results for the profiles and correlation functions (All)-(A13) are in full 



agreement with the exact results of Refs. [11, H3] for large systems, when 
terms of order (a//3) have been neglected. 

In close parallel to section 5.2, we may also calculate the spatial and tem- 
poral correlation functions in the bulk of the low density phase. By setting 
the microscopic order parameter tjTj+i = in (A1)-(A2), the microdynamic 
equation for bulk sites (i <C L — Xr) reduces to 

Ti(t + l)=Ti- 2 (t) ; T l+1 (t) = 0, (A.14) 

where i is even. The correlation function in the NESS with i and R even 
are then, 

< T i+R (t)n(0) > = < n +R _2t(o)Ti(o) > 

= < T i+R (t) >< n >= a 2 , (A.15) 
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holding for < {i,j} « L — Xr with j = i + R — t and R 7^ 2£. For % and/or 
R odd, the correlation function vanish, as < r >= for odd sites. We observe 
that the occupations between two even sites in the low density phase are 
uncorrelated, as a consequence of the sublattice - parallel dynamics. In the 
corresponding case of fully parallel dynamics, the occupations are correlated, 



as is shown in (|43| ) and (45). The absence of correlation is understandable 
here, as a particle, attempting to enter the system in the low density phase 
at site 1 is never blocked. 
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Figure Captions 

Figure 1: Flux j as a function of density p at constant removal rate (3 with 
(3 < 0'. For p F < p < pj with p F = (1/(1 + f3) and pj = 1/(1 + (3) the flux 
j = pp remains constant. The triangle bounding the region j < min{p, 1 — p} 
of coexisting phases is called the fundamental diagram in traffic flow prob- 
lems. 

Figure 2: Instantaneous interface width w between coexisting phases as a 
function of time, measured for L = 1000, a = (3 = 0.9. Notice the saw-tooth 
behavior of the width, corresponding to right and/or left jumps of the in- 
terface boundaries. 

Figure 3: Average interface width < w > between coexisting phases, mea- 
sured over an interval of 2 x 10 4 time steps, as a function of injection rate 
a, in a system with L = 1000, and compared with the theoretical prediction 
< w >= 2/(1 — a). As a | 1, the width of the interface changes from mi- 
croscopic to macroscopic. 

Figure 4: (a) Density profiles for {a = 0.50; (3 = 1} and {a = 0.50 and (3 = 
0.51} for a system with L = 1000. The bulk densities are the same (see (16)), 
except in the boundary layer, where the excess density is < TjTj + i >, aver- 
aged over 1.3 x 10 7 time steps, (b) The 3-point correlation < TjTj + iTi + 2 >. 
In both plots the theoretical and simulation results are indistinguishable. 
The latter are time averaged over 1.8 x 10 4 time steps. 

Figure 5: Density profiles in the coexistence region, time averaged over 4000 
time steps. At time t\/t2 the interface is at R is 302/649, and the average 
number of particles < N > is 566/450 respectively. 

Figure 6: Histogram of travel times for a system of L = 1000 (a) in the low 
density phase with {a = 0.50; (3 = 0.51} and (b) in the high density phase 
with {a = 0.51; p = 0.50}. 

Figure 7: Average travel time < T > as a function of a = (3 for coexisting 
phases in a system with L = 1000 sites. 



Figure 8: 'Equation of state' for the nuxj(p) from (|63|) and (p4| ) (dashed line) 
compared with simulation results for L = 1000 sites at (3 = 0.2(pp ~ 0.17) 
and (3 = 0.5(pf — 0.33). The smooth crossover at p = pp is derived from 
the profile ( |6l| ) of the blockage region. 
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Figure 9: Correlation function < TjTj+iTi + 2 > as a function of i. The points 
are simulation results; the curve represents the relation (pq), with parame- 
ters L = 1000, (3 = 0.5 and p = 0.3. 



Table Captions 

Table 1: Equation of state j(p,f3). 





L= 100 


(5 = 0.25 


PF = 1/5 


N 


P 


j(sim) 


j(theor) 


10 
15 
18 


0.10 
0,15 
0.18 


0.0955 
0.1401 
0.1638 


0.0958 
0.1384 
0.1599 




L= 100 


/3 = 0.5 


PF = 1/3 


20 
30 


0.2 
0.3 


0.1964 

0.2882 


0.1962 

0.2857 




L = 1000 


13 = 0.5 


PF = 1/3 


300 
330 


0.30 
0.33 


0.2981 
0.3243 


0.2980 
0.3230 




L= 100 


(3 = 0.75 


PF = 3/7 


30 
40 


0.30 
0.40 


0.2971 
0.3914 


0.2977 
0.3910 




L = 1000 


= 0.1 


p F = 1/11 


60 

85 


0.60 
0.85 


0.0587 
0.0806 


0.0586 
0.0797 



39 



U.b 
0.4 

0.2 


- 


1 i ' i ' i 


- 


/ /8' \ 
)8 


- 







,i,i, 





0.0 Pp 0.2 



0.4 0.6 

density 



0.8 Pj 1.0 



Figure 1: Flux j as a function of density p at constant removal rate (3 with 
/3 < /3'. For pf < P< PJ with p F = /3/(l+/3) and pj = 1/(1+/?) the flux j = 
/9p remains constant. The triangle bounding the region j < min{p, 1 — p} of 
coexisting phases is called the fundamental diagram in traffic flow problems. 
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Figure 2: Instantaneous interface width w as a function of time, measured 
for L = 1000, a = (3 = 0.9. Notice the saw-tooth behavior of the width, 
corresponding to right and/or left jumps of the interface boundaries. 
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Figure 3: Average interface width (w), measured over an interval of 2 x 10 4 
time steps, as a function of injection rate a, in a system L = 1000, and 
compared with the theoretical prediction (w) = 2/(1 — a). As a | 1, the 
width of the interfacechanges from microscopic to macroscopic. 
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Figure 4: (a) Density profiles for {a = 0.50; (3 = 1} and {a = 0.50 and 
f3 = 0.51} for a system with L = 1000. The bulk densities are the same (see 
(16)), except in the boundary layer, where the excess density is (tjTj+i), 
averaged over 1.3 x 10 7 time steps 
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....(b) The 3-point correlation (tjTj+iTj-)^). I n both plots the theoretical and 
simulation results are indistinguishable. The latter are time averaged over 
1.8 x 10 4 time steps. 
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Figure 5: Density profiles in the coexistence region, time averaged over 4000 
time steps. At time ii/^2 the interface is at R is 302/649, and the average 
number of particles (N) is 566/450 respectively. 
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Figure 6: Histogram of travel times for a system of L = 1000 (a) in the low 
density phase with {a = 0.50; (3 = 0.51}.... 
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.and (b) in the high density phase with {a = 0.51; (3 = 0.50}. 
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Figure 7: Average travel time (T) as a function of a 
phases in a system with L = 1000 sites. 
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Figure 8: 'Equation of state' for the flux j(p) from (63) and (64) (dashed 
line) compared with simulation results for L = 1000 sites. The smooth 
crossover at the p— values pp and pj is derived from the profile (61) of the 
blockage region. 



0.30 



£ 0.24 
o 



u 

I 0.18 



c 

o 



o 

CD 



0.12 



8 0.06 



0.00 
970 




975 980 985 990 995 1000 

site 



Figure 9: Correlation function {TiTi+iTi+2) as a function of i. The points are 
simulation results; the curve represents the relation (66), with parameters 
L = 1000, 13 = 0.5 and p = 0.3. 



